Quantum Fluctuations and Resistivity of Thin Superconducting Wires 
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We present a microscopic study of the quantum fluctuations of the superconducting order pa- 
rameter in thin homogeneous superconducting wires at all temperatures below Tc- The rate of 
quantum phase slip processes determines the resistance R(T) of the wire, which is observable in 
very thin wires, even at low temperature. Furthermore, we predict a new low temperature metallic 
phase below a critical wire-thickness in the lOnm range, in which quantum phase slips proliferate. 



The role of fluctuations for the superconducting tran- 
sition in reduced dimension is well known. Above the 
critical temperature Tc, such fluctuations yield an en- 
hanced conductivity M. Below Tc thermally activated 
phase slips (TAPS) result in a finite resistance of one- 
dimensional (ID) superconductors 0]. Close to Tc the 
experimental results || fully confirm the theoretical pre- 
dictions |^]. However, as the temperature is lowered the 
number of TAPS decreases exponentially and no mea- 
surable resistance is predicted by the theory || at T not 
very close to Tc- Nevertheless, recent experiments by 
Giordano clearly demonstrate a notable resistivity of 
thin superconducting wires far below Tc- More recently, 
similar experimental results on (quasi-)lD systems have 
been also reported by others 

As these experimental findings cannot be understood 
within the TAPS theory it is tempting to attribute 
them to quantum fluctuations which generate quantum 
phase slips (QPS) in a ID superconducting wire. A first 
estimate for the QPS rate oc exp(— Sqps), however, 
leads to the disappointing conclusion that Sqps roughly 
equals the number of transverse channels Nch = k F S 



in the wire (S 



is the cross section of the wire), 



which is very large even for the thinnest wires used 
in the experiments B (e.g. for tq ~ 10~ 6 cm we 
have Sqps ~ 10 2 -j- 10 3 ), and therefore QPS effects 
should be strongly suppressed. This estimate is obtained 
from the formula Sqps ~ Uqps/uq,, with energy bar- 
rier Uqps and attempt frequency uj a ~ A. Assum- 
ing Uqps to be the condensation energy A(0)A 2 /2 in 
a volume during a time A -1 , one obtains Sqps ~ 
£oSN A 2 /2A ~ hpS/Air 2 ~ N C h- A similar estimate has 
been obtained using a phcnomenological time dependent 
Ginzburg-Landau (TDGL) free energy with second order 
time derivatives [Mpl. Furthermore, recently Duan M ar- 
gued that the electromagnetic field yields an additional 
suppression of the QPS rate by the factor exp(— 1/a), 
a = 1/137 is the fine structure constant. Even further 
suppression of this rate — similar to the case of Joseph- 
son junctions M — can be expected due to dissipative 



currents in the QPS core. In contrast, the magnitude of 
the resistance for the thinnest wires measured in Rcf. Q] 
yields Sqps ~ 10 with the QPS rate by orders of magni- 
tude larger than that derived from the above estimates. 

In this Letter we argue that the above estimate needs 
qualitative improvement. First, the estimate for the po- 
tential barrier should be corrected: as the typical electron 
mean free path in the wires [Q is very small I < 10 nm, 
one should rather take £ ~ V^o as the typical QPS size. 
This decreases the estimated value of Sqps by one order 
of magnitude. Second, we show below that the role of the 
electromagnetic field for thin wires was overestimated in 
Ref. § (roughly by a factor r /A L - 10" 1 ^ 10" 2 , X L 
is the London length of a bulk superconductor). Third, 
the dissipative currents turn out not to have a strong 
impact on the QPS rate, especially in the limit of low 
T. Also from a general point of view, TDGL-based theo- 
ries (?]|| are insufficient at T not very close to Tc and fail 
to give qualitatively correct results. A microscopic theory 
of QPS is needed that accounts for non-equilibrium, dis- 
sipative and electromagnetic effects during a QPS event. 

This theory is reported upon below. Taking into ac- 
count the above effects we determine both the character- 
istic size xq and time scale To of a QPS core (in general 
these depend on dissipation and electromagnetic charac- 
teristics of the wire and do not coincide with £ and 1/A) 
and derive a microscopic expression for the QPS rate. 
At T — we predict a new Metal-Superconductor {M S) 
phase transition governed by the strength of electromag- 
netic interactions between different QPS. We evaluate the 
effective resistance of a ID superconducting wire, demon- 
strate that in sufficiently thin wires QPS's are observable 
at all T, and determine the crossover temperature be- 
tween the regimes governed by TAPS and QPS. 

The model. Our calculation is based on the effective 
action approach for a BCS superconductor jl(|. The 
starting point is the partition function Z expressed as 
an imaginary time path-integral over the electronic fields 
ip and the gauge fields V, A, with Euclidean action 
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S = / d A r / dT 



-giplipli/>li>l + ieVrii + [E 2 + B 2 ]/8tt 



Here f3 = 1/T, £(V) = -V 2 /2m - /x describes a sin- 
gle conduction band, 5 is the BCS coupling constant, 
eni denotes the background charge density of the ions, 
and units in which h = kp = 1 are used. A Hubbard- 
Stratonovich transformation introduces the energy gap 
A as an order parameter and the electronic degrees of 
freedom can be integrated out. What remains is an ex- 
pression for the partition function in terms of an effec- 
tive action for A, V and A, with a saddle-point solution 
A = Abcs and V = A = 



P 

(Ft I dr 





|A| 2 E 2 + B 2 



8?r 



TrlnG" 1 , 



O : - ( d T + -{V,v s } ) 1 + Aoar 



HIV, . , . 

£(V) + — ^ - ) a 3 , 



where the superfluid velocity v s = (l/2m)[V<p — 2eA/c], 
the chemical potential for Cooper pairs <& = V — <p/2e, 
and A = Aoe I¥3 have been introduced. 

Effective action. The effective theory is constructed by 
expanding up to second order around the saddle point in 
$ and v s to obtain the electronic polarization terms. Us- 
ing the Ward identity from gauge invariance, the result 
can be written as the sum of terms, related to the normal 
and superfluid densities n n and n s , that describe normal 
and superconducting "screening" respectively |Tl|,p^| 



|w^|>A ' 



dx-^—E 2 

2\U U \ 



s 



dxdT (^ v 2 + e 2 N [^V 2 + ^$ 2 ' 
V 2 n n 



where use was made of the one-dimensional nature r$ < £ 
of the problem, and a is the conductivity of the wire. 
Transverse screening is irrelevant if the London penetra- 
tion depth > r and we retain only one component of 
the vector potential 

A phase-slip event in imaginary time involves a sup- 
pression of the order parameter in the phase slip core 
and a winding of the superconducting phase around this 
core We now separate the total QPS action Sqps 
into a core part Score around the phase slip center for 
which the condensation energy and dissipation by nor- 
mal currents are important (scales x < xq, t < To), and 
a hydrodynamic part outside the core S ou t which depends 
on the hydrodynamics of the electromagnetic fields 



dxdr 



'C + C 



V 2 



c 



1 



2ic 2 



-.A 2 



2e 2 L 



where the kinetic inductance L = raj (e 2 n s S) and kinetic 
capacitance C = Se 2 Non s /n have been introduced, as 
well as C = Se 2 Non n /n which we will drop from now 
on in the limit n s 3> n n at low T. The geometry and 
screening by dielectrics outside the wire are accounted 
for by the capacitance per length C = e r [2 ln(a;o/ro)] _1 
and the inductance times length L = 2h\(xQ / Tq) / c 2 that 
replace the E 2 + B 2 -term. Here c is the velocity of light 
and e r the dielectric constant of the substrate. 

Single QPS. Outside the phase slip core, the equations 
of motion for V, A, and tp are solved by the saddle point 



(p = arg(.x + zc t) 
1 d T ip 



c- 



C~ 



V 



1 + C/C 2e 



A = 



L 

c 



L 



1 + L/L 2e 



The space— time anisotropy is determined by the plas- 
mon velocity cq, rather than by vp. For generic pa- 
rameters the velocity cq reduces to the velocity of the 
Mooij-Schon mode, which has dispersion 0J 2 = c\ [S k 2 
with c 2 MS = Suj 2 p1 /AttC @, where Upi is the 3D plasma 
frequency. The corresponding saddle point action is 



S out = iihi[ mm ( c o/^ X)/ max(c T , x )} 



(1) 



with fi = 7r/[4e 2 c (£ + L)]. 

The contribution from the core part is estimated to be 



q 



,.„- c = ^-Ar x A + - 



E 
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^-\E{^)\ 2 . (2) 



The first part is the condensation energy that is lost in- 
side the core and the second part defines the energy of 
dissipative currents inside the core during a phase slip 
event. Inserting the saddle point solution, and mini- 
mizing the action with respect to Xq and r , we find 
Xq w c t ~ (CTC /2e 2 A f Ag) 1 / 3 and S* ore as three times 
the condensation energy in (§). In the Drude limit 
a = 2e 2 N VFl/3, we obtain x « (£ 2 cq/A ) 1/3 and 



S* orc = 6A 5(Ao0 4/3 /co /3 , 6-1- 



(3) 



QPS Interactions. The next step is to consider a gas 
of QPS's in a superconducting wire. We also assume 
that an applied current I (much smaller than the de- 
pairing current) is flowing through the wire. Substitut- 
ing the saddle point solution tp = t^{ x — x %i t — T 
into the action and keeping track of the additional term 
fdrj dx(I/2e)d x cp @, we find 



S n 



1ST 



v ^w ■ln^w^v 



\Xq 



(4) 



Here pij = (c^(r, — r,) 2 + (pa — xj) 2 ) 1 / 2 defines the dis- 
tance between the i-th and j-th QPS in the (x, t) plane, 
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Ui = +1 (-1) are the QPS (anti-QPS) "charges", and 
$o = hc/2e is the flux quantum. Only neutral QPS con- 
figurations with u to t = v i = ( an d hence n even) 
contribute to the partition function. This is a conse- 
quence of the boundary condition ip(x,r) — ip(x,T + (3) 
in the path integral for the partition function JToj ] . 

Metal- Superconductor phase transition. For 1 = 
Eq. (||) defines the standard model of a 2D gas of log- 
arithmically interacting charges v^. The effective (small) 
fugacity y of these charges (or the QPS rate y/xoTo) is 



y = x T Bexp(-S* OTe ) 



(5) 



where B is the usual fluctuation determinant with zero 
modes excluded. From the Coulomb gas analogy, we con- 
clude that a KTB phase transition 15 for QPS's occurs 
in a superconducting wire at fx = ii* = 2 + Any « 2: 
for \x < ll* the density of free QPS in the wire (and 
therefore its resistance) always remains finite, whereas 
for (i > ll* QPS's and anti-QPS's (AQPS) are bound in 
pairs and the linear resistance of a superconducting wire 
is strongly suppressed and T-dependent. We arrive at an 
important conclusion: at T = a ID superconducting 
wire is in a superconducting state, with vanishing lin- 
ear resistance, provided the electromagnetic interaction 
between phase slips is sufficiently strong, i.e. ll > /i*. 

The above analysis is valid for sufficiently long wires. 
For typical experimental parameters, however, X < cq/3 
(or even X <C cq/3), and the finite wire size needs to be ac- 
counted for. To this end, we first apply the standard 2D 
scaling analysis pl| difx = — 47r 2 yu 2 y 2 and diy = (2 — fx)y, 
where ll and y depend on the scaling parameter I. Solv- 
ing these equations up to I = lx = hx(X/xo) we obtain 
the renormalized values p = fx(lx) and y = y(lx)- For 
larger scales / > lx only the time coordinate matters 
and the problem reduces to a that of a ID Coulomb gas 
with logarithmic interaction. Therefore, (for y <C 1) fur- 
ther scaling for I > lx is defined by jl6|,|l0| d^x = and 
diy = (1 — fx)y. For p > 1 the fugacity scales down 
to zero, which again corresponds to a superconducting 
phase, whereas for p < 1 it increases indicating a resis- 
tive phase in complete analogy to a single Josephson junc- 
tion with ohmic dissipation. Thus, our above conclusion 
about the presence of a MS phase transition at T = 
remains valid also for relatively short wires, although the 
phase transition point is given by the somewhat different 
condition p = 1. In practice, both conditions fx > ll* and 
p = 1 are realized in wires with diameter 2ro ~ 10 -j- 20 
nm, see also the discussion below. 

Wire resistance at low T. At any nonzero T the wire 
has a nonzero resistance R(T, I) even in the "ordered" 
phase [i > li* (or p > 1). In order to evaluate R(T) in 
this phase for a long wire we proceed perturbatively and 
first calculate the free energy correction SF due to one 
bound QPS-AQPS pair. [See Ref. @ (Chapter 5.3) for 
a discussion of a similar procedure.] The one QPS-AQPS 
pair contribution SF to the free energy of the wire is 



5F = 



X V 2 f f X ^ e (<S>oIr/c)~2^n[p(r.x)/x Q ] ^ (g) 

xqt J Tq t J Xq x 



where p = (cqT 2 + a; 2 ) 1 / 2 . It is convenient to first in- 
tegrate over the spatial coordinate x and take the wire 
length X — > oo. For nonzero / the expression in Eq. (|^) 
is formally divergent for (3 — > oo and acquires an imag- 
inary part Im SF after analytic continuation of the in- 
tegral over the temporal coordinate r jl7|,[l0]]. This in- 
dicates a QPS-induced instability of the superconduct- 
ing state of the wire: the state with a zero phase dif- 
ference Sip(X) = tp(X) — (p(0) = decays into a lower 
energy state with Stp(X) = 2tt. The corresponding de- 
cay rate is T2t, = 2lmSF. The rate for the opposite 
transition (which is nonzero at nonzero T) is defined 
analogously with / — > —I. The average voltage drop 
V = (^ /c)[T 2 Tr(I) - r 27r (— I)] across the wire is 



V 



§vXy 2 v^r(M 



ctqXq T(jj,)T(2fj, 



^— sinh ( — 
1) \2cT) 

-i 2^—2 



)' 3 


27TT 







(7) 



r(x) is the Euler gamma-function. For the wire re- 
sistance R(T,I) = V/I this yields R oc T 2 ^" 3 and 
R oc pt 1 -^ for T > $ -f and T < <& I respectively. For 
thick wires with (j, > Li*, we expect a strong temperature 
dependence of the resistivity. For thinner wires the tem- 
perature dependence of the resistivity becomes linear at 
the transition to the disordered phase in which our anal- 
ysis is not valid. At T ®oI/c we expect a strongly 
nonlinear I-V characteristic V ~ I a in thick wires, and 
a universal a(/i*) = 2 in thin wires at the transition into 
the resistive state with V ~ /, i.e. a = 1. Note that in 
contrast to the KTB transition in 2D superconducting 
films, the jump is from a = 2 to 1, instead of a = 3 to 1. 

For a short wire X < cq/T we again proceed in two 
steps. A 2D scaling analysis yields the "global" param- 
eters y, p, and the microscopic cut-off tq = tqX/xq. 
In analogy with the resistively shunted Josephson junc- 
tion [ fl0| , the voltage drop from the imaginary part of the 
free energy reads 



V 



— ; — sinh 

r(2/x)cr \2cT 



|r(A- 



27rfo 



i^nl 



2ttcT [ P 



2/i-l 



giving R oc T 2 ^ 2 and R oc I 2 ^- 2 respectively at high 
and low T. This result is valid for p > 1 and also for 
smaller p at not very small T |icj]. At T — > in the 
metallic phase the resistance flattens off and becomes JTy] 



R = R q /P 



(8) 



where R q = 7r/2e 2 ~ 6.5 kf2 is the quantum resistance. 
Note that in this limit the size dependence of the resis- 
tance enters only through the renormalized value p. 
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Crossover temperature. The cross-over between the 
TAPS and the QPS regime occurs at a temperature T* 
that can be estimated by comparing the corresponding 
resistivities Rtaps and Rqps- The standard result of 
Ref. g is Rtaps = f3n(h/4e 2 )exp{-f3AF), w ith A F = 
V2HfSZ/3TT, attempt frequency fi = (X/Ov^STt^ 1 , 
and relaxation time r^ 1 = 8(Tc — T)/ir. For thin wires 
with /i < S* OIe the value T* which follows from a com- 
parison of the exponents (3AF and 25* ore is 

T* = £L « (9) 

core 

The pre-exponential factor i? in Eq. ([5]) can be estimated 
by matching the pre-exponential factors of Rqps and 
Rtaps at T = T*. We postpone a more detailed analy- 
sis of the pre-exponent to a forthcoming publication jl2) . 

Discussion. For typical system parameters A:^ 1 ~ 0.2 
nm < Z ~ 7 nm < £ ~ 10 nm < £0 ~ Al ~ 100 nm, we 
find that L and C drop out of the problem and L and 
C determine the physics (unless T ~ T c or e r 3> 1, see 
Ref. p2[ ). We will also take the length X of the wire to 
be smaller than the localization length, so that localiza- 
tion effects do not play a role. Taking ro ~ 10 nm and 
e r = 1, we obtain the velocity cq/c — cms /c ~ (ro/6Az,), 

/i = {n^;/8a)(r /\ L ) « 50(r /A L ) , 

and 2S* ore < 10. This estimate demonstrates that — in 
contrast to previous studies j|] — quantum fluctuations 
in thin superconducting wires are not negligibly small 
and can be well observed in experiment. Furthermore, 
our estimate for the classical-to-quantum crossover tem- 
perature Eq. (g) yields T* ~ 10A(T*), i.e. for thin wires 
one expects this crossover to happen quite close to Tc- 
These features are in good agreement with the experi- 
mental findings . 

For the quoted parameters, we predict the supercon- 
ductor to metal transition at a wire thickness of order 
ro = Ai/25 « 5 -j- 10 nm. This prediction agrees with the 
results of Giordano, who finds that wires with ro ~ 8 nm 
have a resistivity that saturates at a measurable level at 
low T, whereas the resistivity of thicker wires ro ^ 13 nm 
decreases with temperature even at the lowest tempera- 
tures {§ . Also the saturation value R ~ 10 ~ 20 kfl (g) 
predicted here is consistent with that measured in j| . 

Note that superconductivity in wires with radius ro ~5 
nm is not destroyed by finite size and level spacing ef- 
fects; particles of radius down to 1-3 nm turn supercon- 
ducting p8l. Finally, the MS phase transition discussed 
here is in many respects different from that in granular 
wires |l|] . In the latter case the role of disorder and the 
wire thickness is irrelevant, and it is the onsite Coulomb 
interaction that drives the transition into an insulating 
phase. For homogeneous wires, in contrast, the transition 
is into a metallic state. 



In conclusion, we have studied QPS's starting from mi- 
croscopic theory and find a measurable resistivity in su- 
perconducting ultrathin wires at temperatures T <C Tc, 
as well as a new superconductor to metal phase transition 
as a function of the wire thickness. 
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